%求解控制序列为p时候的相关参数
%t_f = x(1)、c_1 = x(2)、c_2 = x(3)
%编写一个计算 F 的函数，即 x 处方程的值
            %中间变量t1时刻的速度与位移
            F = @(x)[12-0.5*x(2)*x(1)*x(1)+x(3)*x(1)-12*(1-0.15*sin(0.2*x(1)));
                30-0.15*x(2)*x(1)*x(1)*x(1)+0.5*x(3)*x(1)*x(1)+12*x(1)-(2*10-2*1.8*12*(1 - 0.15 * sin(0.2 * x(1)))+(100 + 12 * x(1) + 10 * cos(0.2 * x(1)) - 10))];
            
%创建初始点 x0 
x0 = [0;0;0];
%设置选项以返回迭代显示。
options = optimoptions('fsolve','Display','iter');